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ABSTRACT 



A problem of routing earth resource survey 
aircraft, proposed by NASA, is formulated as a 
traveling salesman problem, in which the salesman 
(aircraft) has a range constraint. A heuristic 
algorithm is presented, which seeks a minimal length 
set of subtours through n cities. The aircraft begins 
a subtour at the base location, visits a subset of the 
n cities and returns when the range constraint 
prevents a visit to another city. Additional subtours 
are created until all cities are visited. 

The algorithm is programmed in FORTRAN for use on 
digital computers. The IBM 360/67 computer at the 
Naval Postgraduate School was used to find solutions 
to three operational problems of size seven, eighteen 
and twenty-five cities. Computation times for each 
problem was under 20 seconds and the solutions were 
significiantly better than feasible solutions 
calculated without the use of the algorithm. 



4 



TABLE OF CONTENTS 



I. I NTBODUCTIO N 7 

II. PROBLEM FOBMU LATION 10 

III. SOLUTION TECHNIQUES 15 

IV. THE ALGORITHM 19 

V. TEST EROBLEMS 22 

VI. ADDITIONAL APPLICATIONS 32 

VII. CONCLUSIONS 34 

Appendix A: COMPUTER PROGRAM 36 

LIST OF REFERENCES 37 

INITIAL DISTRIBUTION LIST 48 



5 



ACKNOWLEDGEMENT 



It is with pleasure that I wish to acknowledge Roger D. 
Arno, Research Scientist at the National Aeronautics and 
Space Administration* s Ames Research Center, Moffett Field, 
California for his assistance in formulating this problem 
and for providing some of the data used in the examples. 



6 



I. 



INTRODUCTION 



The National Aeronautics and Space Administration is 
fostering many technology utilization programs to explore 
the application of its aerospace technology' to other fields. 
Suggesting how our finite natural resources can be used most 
effectively is one of its more important considerations. 
NASA has been working to sophisticate the remote sensing of 
the earth from a variety of aerospace platforms. The 
information available from such a program can provide 
accurate and timely data to assist in the management of the 
earth’s natural and cultivated resources. NASA has 
demonstrated the usefulness of remote sensing to a wide 

variety of users with emphasis on the disciplines of 
geology, agriculture, forestry, oceanology, meteorology and 
hydrology. For example, an agriculturist can be provided 
with identification and area measurements of major crop 
types for use in forecasting potential crop yield. In 

addition, remote sensing can sometimes detect agricultural 
diseases and pest attacks before they can be detected from 
the ground. More detailed discussions of the type of 

instruments used and examples of the data available from 
remote sensing can be found in references 4, 6, 7, 13, 15, 
and 16. 

In recent years land has emerged as the one resource 
"basic to both the good life and the goods of life" [ref 8]. 
There has been emphasis on efficient land use planning, 
particularly in metropolitan areas. Many experts are 

convinced that a well conceived program of remote sensing of 
metropolitan areas will enhance the quality of urban design. 
This program could provide, for a given area of land, its 
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access 



biological productivity, soil limitations, topology, 
to transportation and utilities, existing vegetation and the 
use of adjacent land areas. This data would be a very 
valuable aid in efficient urban planning. 

NASA has been involved in the development of 

instrumentation to perform the variety of sensing required 
and demonstration of the data available to potential users. 
In addition, NASA has been asked by potential users to 

estimate total program cost and to recommend base locations 
for a special fleet of aircraft to support a program of 
remote sensing of urban centers. Assuming an initial cost 
for the program set up and instrument procurement and a 
fixed cost associated with each sensing mission, NASA has 
been seeking methods of analyzing the reduction of operating 
costs through efficient vehicle routing and effective base 
locations. 

The purpose of this paper is to develop an algorithm 
which will, given a base location and a set of urban areas 
to be sensed, provide an optimal or near optimal routing for 
the number of tours required to visit each of the areas. It 

is assumed that the operating cost will be primarily a 

factor of distance traveled plus a fixed charge for each 
tour required to visit the set of urban areas. Using the 
fewest flights (subtours) to visit each area and optimally 
routing the aircraft to insure the minimum distance traveled 
will yield the minimum operational cost. 

This problem was formulated as a multiple traveling 
salesman problem with the addditional constraint of aircraft 
range. The aircraft or sensing platform may only be capable 
of visiting a subset of the set of urban centers before it 
must return to its base (completing a subtour) to refuel. 
Thus the problem is to minimize the sum of the distances of 
the subtours. The algorithm presented by Raymond [ref 14] 
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was modified to consider this range constraint. The 
modified algorithm was programmed in FORTRAN and will be 
discussed along with operational results and recommendations 
for additional uses. 
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II. PROB LEM FORMULAT ION 



The traveling salesman problem (TSP) is one of the 
classic challenges in operations research (OR) . The problem 
is simply stated: a salesman starting from his home wishes 
to visit a set of cities once and only once and return home. 
In what order should the cities be visited to minimize the 
total distance traveled? Mathematical formulation of the 
problem and a discussion of the constraints can be found in 
Wagner [ref 17] Chapter 13. 

OR literature has included applications of this model to 
a variety of practical problems. Industrial organizations 
have used the model to improve the routing of delivery 
vehicles, an important consideration in light of increased 
oil prices and rising driver salaries. The "Chinese postman 
problem" is a similiar problem applied to postal service 
introduced by the Chinese mathematician Kwan [ref 10]. 
Municipal governments have found the model valuable for 
routing garbage collection, street cleaning, meter reading 
and snow removal vehicles. Other applications include dairy 
and newspaper delivery. Further discussion of applications 
can be found in Bradley's "Survey of Deterministic Networks" 
[ref 3], 

The key difference in the classic statement of the TSP 
and the remote sensing problem is the addition of a range 
constraint. The salesman of the TSP may refuel during his 
trip without altering the optimal route. This is not 
feasible for the aircraft performing remote sensing. There 
exists a point during the trip when it must return to its 
base to refuel and then continue the mission. Therefore, 
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the problem statement becomes: an aircraft starting from its 
base must visit a set of urban areas once and only once and 
return to -its base. The aircraft is subject to a range 
limitation for any given flight; that is, it may only be 
capable of visiting a subset of the set of cities be-fore it 
must return to its ba-Se to refuel. Each time the aircraft 
returns to its base a subtour is completed and additional 
flights (subtours) are created as necessary to accomplish a 
visit to each city. The objective is tc minimize total 
distance traveled as the sum of the length of each subtour. 

A subtour, t, can be represented as a set of ordered 
city pairs, e.g. 

t = [(i#i)#(i,i) • • • (i )#(i t i ) 1 

12 z 3 n-1 n n 1 

where i is always the base. The mathematical formulation 
1 

of the problem is to find the set of subtours to: 



minimize 



E 



k=1 



2 . 1 



subject to 1 = 

k 



= ’y* a. . ]♦ «. . 

j = 1 K 



2 . 2 



1 < r k=-1 ,2,. . . ,s 

k 



2.3 



1 < s < nc 



2.4 



each city in exactly one subtour 



11 



where s = number of subtours 

1 = length of subtour k 

k 

n = number of cities in subtour k 

k 



d = distance from city i to city j 

i j 

r = range of aircraft 

nc = number of cities to be visited 

Note that if s=1, eg. 2.2 replaces eg. 2.1 and the problem 
is a standard TSP. The upper bound of eg. 2.4 represents 
the worse case situation that each city must be visited by a 
seperate subtour. Eg. 2.3 represents the range constraint. 

The distance matrix (D ) represents the euclidean 

ij 

distance from city i to city j. During actual sensing 

missions the aircraft must fly several passes over a city to 

completely image the area. This distance (loiter distance) 

depends on three factors: the size of the area, the percent 

side overlap used for film images and the turning radius of 

the aircraft. Each city will have a unigue loiter distance 

(Id ) that must be considered when the subtour length is 
i 

calculated. When the problem is formulated the value Id is 

i 

added to each element in column i. Thus the new distance 

matrix is D' = D + Id for all i. The new matrix is 
ij ij i 

used to compute subtour lengths. 



This problem appears similiar to two existing expansions 
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of the TSP; the multisalesman problem (MTSP) and the vehicle 
dispatching problem. In the MTSP, m salesmen starting from 
a base city must visit n cities and return. Each city is 
visited by one and only one of the m salemen. The objective 
is to route each of the m salesmen such that the total 
distance traveled is minimized. Bellmore and Hong [ref 2] 
provide a detailed discussion and formulation of this 
problem. The objective of the vehicle dispatching problem is 
to obtain delivery routes from a depot to a set of demand 
points for a fleet of vehicles such that the total distance 
traveled by the entire fleet is minimized as with the m 
traveling salesmen. The problem is enhanced by finite 
vehicle capacities, route time constraints and varying 
quantity demanded at each demand point. A complete 
description cf this type application can be found in Golden 
[ref 12]. 

There are important differences between the remote 
sensing problem and the applications mentioned above. The 
MTSP has no upper bound on the number of cities any one of 
the m salesmen can visit. The aircraft range constraint 
requires consideration of a upper bound on the number of 
cities that can be visited during any given flight. The 
vehicle dispatching problem is concerned with routing each 
vehicle in the fleet. Although potential users will 
presumably have a "fleet" of aircraft to accomplish its 
sensing missions, not every aircraft will be scheduled to 
sense a subset of the areas to be visited. For example, 
given a fleet of three aircraft and the determination that 
two subtours are required to visit the set of cities, only 
two aircraft (or one aircraft making two flights) will be 
scheduled. This decision assumes the fixed cost charged 
each time an aircraft is launched is sufficiently high to 
make it more desireable to create the fewest number of 
flights. 
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The remote sensing problem, therefore, is not a simple 
restatement of an existing TSP. The addition of a range 
constraint requires special considerations. The solution 
technique must carefully create only the number of subtours 
required and route each subtour so that the total distance 
traveled is minimized. 
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III. SOLUT ION ALGORITH MS 



One of the most intriguing aspects of the TSP is, that 
although simply stated, it remains unsolved in a closed 
form. The computational difficulty with known solution 
tecnnigues is the exponential increasing solution time as 
the size of the problem increases. There are two types of 
solution algorithms: heuristic and exact. The latter 
yields the optimal tour. The difficulty with this type 
algorithm is the large computer storage and time 
requirements which are related to the large number of 
possible solutions. For example, given a symmetric TSP 
there are (n-1) !/2 possible solutions (n=number of nodes). A 
relatively small problem of ten nodes yields 181,440 
solutions. Fortunately, many applied or real world problems 
do not demand the optimal solution. Powerful heuristic 
algorithms have been developed which will yield near optimal 
solutions and provide a substantial savings in computer 
storage, time and money. An excellent survey of existing 
solutions is given by Bellmore and Nemhauser [ref 2]. 

The Bellmore and Nemhauser article classifies the 
algorithms by solution generation method. There are three 
fundamentally different ways of generating solutions: 
tour-to-tour improvement, tour building, and subtour 
elimination. The remote sensing problem requires tour 
length feasibility be maintained as the tour • is formed. 
That is, if visiting the next city (node) will result in 
violating the range constraint, a new subtour is started. 
Consequently, a tour building solution was considered the 
best technique to adapt to this problem. 
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Both exact and heuristic algorithms were studied for 
applicability. A literature survey provided a number of 
both types of algorithms. An exact solution considered was 
presented by Christofides and Eilon [ref 5]. Their paper 
discusses three solution methods to the vehicle-dispatching 
problem. The branch and bound approach, a modification of 
the algorithm published by Little et. al. [ref 11], would 
reach an exact solution to this problem provided the number 
of subtours required was known a priori. If s (the number 
of subtours) was known, the base location would be replaced 
by s artificial bases as explained in ref 5. The algorithm 
would progress using the branch and bound method to insert 
nodes into one of the tours. Before branching to a new node, 
it would be necessary to check that the total distance 
accumulated on each subtour did not exceed the range. The 
key to this algorithm is knowing the number of subtours a 
priori. 

As noted in the problem formulation section, there 
exists an upper bound on the number of subtours. Therefore, 
it is feasible to select an s known, a priori to be 
sufficiently large, and apply this branch and bound 
algorithm. The solution would be recorded, the number of 
subtours decreased by one and a new soluton computed. This 
procedure would continue until there are insufficient 
subtours to complete a visit to each city and the last 
recorded solution would be the optimal. The success of this 
procedure would be very dependent upon the initial selection 
of s. A poor selection could involve many iterations and 
would require considerable user interaction and computer 
time. This would be particularly undesireable for problems 
that require a solution in a reasonably quick time. 
(Problems of this type are discussed in Section VI) . 
Therefore, application of this algorithm was not considered 
profitable. 
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Limiting the study of heuristic algorithms to those 
using the tcur building generation technique, the algorithm 
presented by Raymond [ref 14] was selected to modify and 
apply to this problem. The approach of Raymonds' s algorithm 
is similiar to that of Karp and Thompson [ref 9] and will 
obtain optimal or near optimal solutions to the classical 
TSP . The algorithm successively inserts nodes into the tour 
by using a series of decision rules and provides for an 
exchange procedure to determine if the total distance 
traveled may be reduced by rearranging the order in which 
the nodes are visited. A detailed explanation of the 
algorithm steps is available in ref 14. 

The modified algorithm starts by creating a subtour from 

the base to the most distant city. The same criteria as the 

original algorithm is used to select the next node for 

insertion. That is, increment calculations (I ) are made 

k 

for each remaining node (k) . This calculation is made by 

inserting the node k between each pair of nodes (p,q) in the 

subtour and computing the incremental increase in subtour 

length, I = d + d - d Each node will have a minimum 

k pk kq pq 

increment calculation (I* ) and a next best increment 

k 

calculation (J* ) . If the number of cities in the tour is 

k 

less than four the maximum I* determines the node k to be 

k 

inserted between the nodes that produced I* . Otherwise, 

k 

the node k for which J* - I* is a maximum is selected and 

k k 

inserted between the nodes that produced I* . Before 

k 

completing the insertion, the algorithm looks ahead to 

determine if the range constraint will be violated. If it 
is not violated the node is inserted and the three link 
exchange procedure described by Raymond is applied. This 
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procedure takes each node in turn out of its present 
subtour, reevaluates its insertion increments in that 
subtour and inserts the node between the nodes that now give 
the minimum tour length increment. If the constraint is 
violated the node is not inserted and an additional subtour 
is formed. 

The nearest city to the base not yet in a subtour is 
used to form the next subtour. The criteria for inserting 
nodes for the multiple tour problem is changed. A node is 
inserted in the subtour that will create the least increase 
in total distance (sum of subtour lengths) . A procedure was 
developed that exchanges nodes among subtours to determine 
if a savings in total distance traveled is possible by 
placing a node in a different subtour. Before any node 
insertions or exchanges among subtours is completed, the 
algorithm looks ahead to insure tour length feasibility is 
maintained. The algorithm progresses iteratively until each 
city is visited. 
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IV. THE ALGORITHM 



The steps of the algorithm: 

1. With the base as one node select the node farthest from 
base (assumed to be j*) and form a tour with links (3,j*) 
and (j*,B). Set the number of subtours s equal to 1 and the 
number of links 1 equal to 2. 



2. Define I* as the minimum tour-length increment that 
k 

results from the insertion of node k in each existing link 

(p,q) . Define nodes p* and q* as the nodes connected by the 

link (p*, q*) that produced I* . In addition define J* as 

k k 

the increment for the node k most nearly equal to I* . Thus 

for node k, I* is the best increment, J* the next best 

k k 

increment and (p*,q*) the link to break if inserting node k. 

Calculate I* and J* for each node that is not currently 
k k 

included in the tour. 



3. If 1 < 4 choose the node m for which I* is a maximum. 

1 k 

If 1 >4 choose the node m for which J* I* is a 

1 k k 

maximum. 



4. Insert the node m into the tour between p* and q* by 
adding links (p*,m) and (m,q*) and deleting link (p*,q*). 
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Calculate the new tour distance d^. If d^ < range, complete 
the insertion and increase the number of links in the tour 
by one. If d^ > range, do not insert node m and go to step 
7 to create a new subtour. 

5. Take each node in turn out of the present tour, 
reevaluate its insertion increments and insert the node 
between the pair of nodes that now gives the minimum tour 
length increment. (In most cases the node is reinserted in 
the link created when the node was taken out of the tour, 
producing no change.) 

6. If all nodes have been added to the tour stop, otherwise 
go to step '2. 

7. Select the node m*, not already in a subtour, closest to 

the base and form a tour with links (B,m*) and (m*,B). Set 

the number of subtours s equal to s+ 1 and the number of 

links 1 equal to 2. 
s 

8. Calculate I* as in step 2 for each node k remaining to 
be added to a subtour. Consider the subtours in the order 

in which they were created. For each subtour select the 

node that produced the minimum I* . Before inserting the 

k 

node two checks are made. First insure that the range 

constraint will not be violated. Second insure that the 
total distance traveled (sum of subtour lengths) will not be 
less if the selected node is inserted in a different 
subtour. Only those nodes which pass these checks are 
inserted. Thus at any given iteration a node need not be 
inserted in each subtour. A given node may have a smaller 
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would force range 



I* in one subtour, but if inserted, 
k 

constraint violation. In this situation place the node in 

the next best subtour. If nodes may not be added to any 
existing subtour without violating the range constraint go 
to step 7. 

9. Repeat step 5 for each subtour to which a node was added. 

10. For each subtour make the increment calculations using 
the nodes from the other subtours. Select the minimum I* 

k 

and check if the total distance traveled may be reduced by 

adding the selected node to the given subtour and deleting 
it from its original subtour. If the distance is reduced 
make the change and repeat step 5. 

11. If all nodes are contained in a subtour stop, otherwise 
go to step 8. 
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V. 



TEST PROBLEMS 



The computer program was used to provide routes from 
bases at Moffett Field, California, Kansas City and Wallops 
Station, Maryland to fifty selected urban centers. Figure 1 
is the grid system used to calculate the distances between 
cities. The cities were placed in one of three sections by 
the following criteria: the straight lines from Moffett 
Field to Kansas City and Wallops Station to Kansas City were 
bisected by a perpendicular, as shown in Figure 1, creating 
three sections. The cities would be imaged from the base in 
their section. Table 1 provides the coordinates and loiter 
distance of each city. The range of the aircraft was 
assumed to be 2900 statute miles. Tables 2-1, 2-2 and 2-3 
are the routing recommendations of the computer program. 

Tables 3-1, 3-2 and 3-3 are the results of a feasible 
solution to the same problem computed at NASA’s Ames 
Research Center. The constraints observed for determining 
this solution were: the aircraft would have a range of 2900 
statute miles, no more than three urban centers would be 
imaged during any subtour and flights from Moffett Field and 
Wallops Station would be held to a reasonable minimum. 

The NASA solution required 29 subtours and a total 
distance of 52,660 statute miles. The algorithm solution 
constructed 15 subtours with a total distance of 36,505 
statute miles traveled. A savings of 16,155 statute miles 
results from the algorithm solution. This large savings is 
due, in part, to the differing criteria used to formulate 
the solutions. However, it does show the usefulness of the 
algorithm and the speed at which the program can compute a 
final solution. 
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X-Y COORDINATE SYSTEM 



TABLE 1 

CITY COORDINATES AND LOITER DISTANCE 



city 


X 


Y 


Loiter 

Distance 


Albany 


2230 


1390 


190 


Allentown 


2190 


1230 


190 


Atlanta 


1950 


700 


290 


Baltimore 


2190 


1150 


340 


Birmingham 


1820 


650 


340 


Boston 


2350 


14 70 


280 


Buffalo 


2030 


13 10 


240 


Chicago 


1660 


1120 


630 


Cincinnati 


1830 


1015 


350 


Cleveland 


1930 


1160 


410 


Columbus 


1900 


1080 


230 


Dallas 


1350 


500 


740 


Dayton 


1820 


1080 


250 


Denver 


910 


900 


450 


Detroit 


1850 


1230 


350 


Greensboro 


2150 


910 


290 


Hartford 


2300 


1350 


150 


Houston 


14 70 


350 


670 


Indianapolis 


1730 


1030 


400 


Kansas City 


1380 


880 


380 


Los Angeles 


220 


560 


670 


Louisville 


1780 


910 


190 


Memphis 


1560 


710 


220 


Miami 


2300 


3 70 


370 


Milwaukee 


1620 


1210 


270 


distances provide 


by Ames 


Research 


Center. 
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CITY COORDINAT 
City 

Minneapolis 
Nashville 
New Orleans 
New York 
Norfolk 
Oklahoma City 
Omaha 

Philadelphia 
Phoenix 
Pittsburgh 
Portland 
Providence 
Rochester 
Sacram en to 
St Louis 
St Petersburg 
Salt Lake City 
San Antonio 
San Bernardino 
San Diego 
San Francisco 
Seattle 
Syracuse 
Toledo 

Washington DC 
Moffett Field 
Kansas City 
Wallops Station 



1 (con'd) 

LOITER DISTANCE 

Loiter 



Y 


Distance 


1260 


380 


800 


270 


4 10 


290 


1300 


730 


2020 


230 


660 


330 


1010 


230 


1200 


480 


4 70 


250 


1150 


430 


1250 


450 


13 70 


160 


1330 


150 


850 


410 


910 


500 


410 


250 


910 


200 


270 


290 


560 


190 


480 


510 


780 


530 


1420 


500 


1350 


310 


1150 


230 


1100 


350 


760 




880 




1 080 





TABLE 

ES AND 

X 

140 0 

1780 

1740 

2300 

2270 

1290 

1300 

2230 

560 

2030 

150 

2360 

2090 

150 

1590 

2140 

580 

1270 

310 

320 

100 

190 

2160 

1830 

220 0 

120 

1380 

2270 
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TABLE 2-1 



FLIGHT MISSIONS FOR URBAN CENTER COVERAGE 
BASE: MOFFETT FIELD 

btour Cities Imaged 

1 Sacramento - Seattle - Portland 

2 Los Angeles - San Diego - San Bernadino 

3 San Francisco - Salt Lake City - Phoenix 

Total Distance = 7772 miles 

Time to Compute Solution 3 1.27 seconds 



Length 

2892 

2184 

2696 
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TABLE 2-2 

FLIGHT MISSIONS FOR URBAN CENTER COVERAGE 
BASE: KANSAS CITY 

btour Cities Imaged 

1 Houston - San Antonio - Oklahoma City 

2 Dallas - Denver 

3 Omaha 

4 Indianapolis - Louisville - St. Louis - Kansas City 

5 Nashville - Birmingham - New Orleans - Memphis 

6 Minneapolis - Milwaukee - Chicago 

Total Distance = 14092 

Time to Compute Solution = 6.97 seconds 



Length 

2878 

2853 

581 

2520 

2746 

2514 
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TABLE 2-3 



FLIGHT MISSIONS FOR URBAN CENTER COVERAGE 
BASE: WALLOPS STATTION 



ubtour 


Cities Imaged 




Length 


1 


St Petersberg - Miami 


- Norfolk 




2642 


2 


Columbus - Cincinnati 


- Dayton - Detroit 


- Buffalo 


2816 


3 


Hartford - Providence 


- Boston - Albany - 


Philadelphia 


2401 


4 


New York - Syracuse - 


Rochester - Allentown - Baltimore 


2607 


5 


Washington/ DC - Pittsburgh - Cleveland - 


Toledo 


2451 


6 


Greensboro - Atlanta 






1724 






Total Distance = 1464 1 

Time to Compute Solution = 17.02 seconds 
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TABLE 3-1 



FLIGHT MISSIONS FOR URBAN CENTER COVERAGE 
BASE: MOFFETT FIELD 



obtour 


Cities Imaged 


Length 


1 


Salt Lake City 


1950 


2 


Phoenix 


2340 


3 


Denver - Omaha 


1990 


4 


Houston 


2020 


5 


San Antonio 


1990 


6 


Dallas 


1990 


7 


Oklahoma City - Kansas City 


1410 


8 


New Orleans - Memphis 


1960 


9 


St Louis - Minneapolis 


1850 


10 


Chicago - Milwaukee 


1800 


1 1 


Miami 


2780 


12 


St Petersburg 


2250 


13 


Birmingham - Atlanta 


2070 


14 


Nashville - Louisville - Indianapolis 


1950 


15 


Cincinnati - Dayton - Columbus 


2100 


16 


Detroit - Toledo 


1880 


17 


Cleveland - Pittsburgh 


2380 


18 


Buffalo - Rochester 


2190 


19 


Syracuse - Albany 


2710 



Total Distance = 39520 miles 
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TABLE 3-2 



FLIGHT MISSIONS FOR URBAN CENTER COVERAGE 
BASE: KANSAS CITY 



S ubtour 


Cities Imaged 


Length 


1 


Seattle 


1950 


2 


Portland 


1450 


3 


San Francisco - Sacramento 


1140 


4 


Los Angeles - San Bernadino 


1510 


5 


San Diego 


1400 



Total Distance = 7450 miles 
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TABLE 3-3 



FLIGHT MISSIONS FOR URBAN CENTER COVERAGE 
BASE: WALLOPS STATION 



Subtour Cities Imaged Length 

1 Boston - Providence - Hartford 1430 

2 New York 1120 

3 Philadelphia - Allentown 1080 

4 Washington, DC - Baltimore 970 

5 Norfolk - Greensboro 1090 



Total Distance = 5690 miles 
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VI 



ADDIT ION AL APPLICA TIONS 



The purpose of this paper was to develop an algorithm 
that will provide optiimal or near optimal routing for a 
range constrained vehicle visiting a set of nodes. The 
algorithm may be used to solve similiar problems. Several 
considerations for use of this algorithm are discussed in 
this section. They are: under-flying satellite coverage, 
isolated disaster coverage and point reconaissance . 

Satellites are often used to image large land masses and 
bodies of water. These orbiting platforms perform the 
sensing mission during periodic passes over the critical 
area. The product of a mission may sometimes be impaired by 
cloud cover or may reveal the need for finer resolution 
imagery. The areas requiring additional coverage are 
assumed dispersed and relatively small. These areas can be 
considered as the cities of the original problem and the 
algorithm used to determine, from a given base, the number 
of flights required, the optimal or near optimal routing and 
the total distance traveled. 

Disasters such as tornados or flooding rivers may 
require government agencies to estimate damage and make 
recommendations for appropriating the relief effort. Remote 
sensing can provide valuable data to assist in the decision 
process. Consider a flooding river; it is likely several 
areas along the flood path need immediate attention. The 
goal is to provide the data for evaluation as rapidly as 
possible. The algorithm can provide routing from a base to 
the areas of concern such that the sensing mission is 
completed in the shortest possible time. It is assumed that 
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least distance routing implies the least time will be 
required to accomplish the mission. 

Several other applications are similiar to the disaster 
coverage and involve routing reconnaissance aircraft. For 
example, when fighting large forest fires, "hot spots" 
develop that need periodic monitoring. Also, in a large 
theater of operations, the military commander often requires 
timely reconnaissance of key objective areas. In both these 
situations the goal is to provide the information to the 
user as rapidly as possible. Therefore, the algorithm can 
be used to route the aircraft to the critical areas and 
provide the information as quickly as possible. 
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VII. CONCLUSIONS 



The algorithm and computer program presented can 
successfully construct an optimal or near optimal solution 
to the range constrained vehicle routing problem formulated 
in this paper. 

Several aspects of remote sensing, considered beyond the 
scope of this. paper, provide a basis for further study in 
this area. The angle of the sun can have an appreciable 
effect on the resolution of the imagery. Sun-angle is a 
function of the time of day and could be considered in the 
computer program by adding a time dimension. Using aircraft 
speed, it would be possible to calculate a time period 
during which the aircraft would image a given city. This 
would provide the user the ability to insure that the 
desired sun-angle occurs during the imaging. 

This paper assumes each aircraft has sufficient sensor 
capacity to image any number of cities in a subtour. It is 
likely there may be a finite capacity for certain sensors. 
For example, assume only a fixed amount of film can be 
placed on the aircraft. The film requied to image all the 
cities in a given subtour may exceed this amount. To 
consider this problem, the aircraft would start with a 
finite amount of film which would be reduced by the amount 
required to image a given city as that city is visited. 
Thus sensor capacity could be added as a constraint and 
considered in the same manner as the range constraint. 

The problem of imaging areas as widely dispersed as the 
fifty urban centers listed in Table 1 requires multiple 
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bases. The example problem assigned cities to bases by 
section. This approach is reasonable, however, it would be 
preferable to permit the computer program to make this 
assignment. One consideration is to construct an artificial 
base from which all flights originate and sufficient dummy 
bases at the same location as the actual bases to permit 
multiple subtours through that base. The distance matrix 
must be constructed such that each subtour starts from the 
artificial base to one of the actual bases and is routed 
back through the actual base to the artificial base. 
Formulation of a fifty city problem with artificial and 
dummy bases could require considerable computer storage. 
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a####*######*#######***#########*##*#*##*#### 



APPENDIX A 



COMPUTER PROGRAM 



VARIABLE DEFINITIONS 

X(I)=X COORDINATE OF THE I-TH CITY 
Y( I ) =Y COORDINATE OF THE I^TH CITY 
Z( I ) =LOI TER DISTANCE OF THE I-TH CITY 
D(I, J)=DISTANCE MATRIX 
RANGE=AIRCRAFT RANGE 
NTOURS=NUMBER OF SUBTOURS CREATED 
LINKS! I )=NUMBER OF LINKS IN SUBTOUR I 
INNODE (It J)=VECTOR OF CITIES IN SUBTOUR I 
NONODE ( I)=VECTOR OF CITIES TO BE VISITED 
LENGTH (II =LENGTH OF SUBTOUR I 
TOTDlS=TOTAL DISTANCE TRAVELED 

LH(I,J)=LEFT HAND ENDPOINT OF SUBTOUR It LINK J 
RH( I t J ) =R I GHT HAND ENDPOINT OF SUBTOUR It LINK J 
NC=NUM BER OF CITIES IN THE PROBLEM 
S(I,J)=IN SUBTDUR I THE SUCCESSOR OF CITY J 
PRED(ItJ)=IN SUBTOUR I THE PREDECESSOR OF CITY J 



STAR(I , J) 
NNSTR( I,J) 


> 




> 




LLSTR ( I , J ) 


> 




SECSTR (I ,J) 


> 


STORAGE MATRICES FOR ALGORITHM 


SNNSTR ( I , J ) 


> 


calculations 


SLLSTR ( I , J) 


> 




THETAt I , J) 


> 





SMS#######*##*#:*# 



C MAIN 



COMMON 0(25,2 5) , S ( 25 , 25 ) , I NNODE ( 25 ,2 5) , NONODE (25) , 

1RH( 25,25) ,LH( 25 , 25 , L I NKS ( 25 ) , THATA ( 25 , 25 ) ,NTOURS,KK, 
HENDtNC ,TOTD IS. LENGTH (25) , RANGE, I STOP, P RED ( 25 , 2 5 ) 
DIMENSION X < 2 5) ,Y(25) ,STAR( 25,25) ,NNSTR ( 25 ,25 ) , 

ILLSTR (25,25) , SECSTR(25 ,25) , SNNS TR ( 25 , 25 ) , SLLSTRt 25,2 5) 
1 ,CRIT(25) , OLD RH( 25,25) , OLDLH ( 25 , 25 ) , Z( 25 ) 

INTEGER S,RH,OLDRH, OLDLH, EXNODE , PRED 
REAL LENGTH 
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INITIALIZE ALGORITHM VARIABLES 



CALL IN I T 
I PR INT = 0 
I F ( I STOP .NE. 



3) GO TO 100 



C 

C 



THE CRITERIA TO FORM THE FIRST SUBTOUR IS TO SELECT THE 
CITY FURTHEST FROM THE BASE AS THE FIRST CITY TO VISIT 



AMAX) GO TO 18 



C 

C 



AMAX=-1 .0 
DO 17 J=1,NC 

IF ( D ( 1 1 J ) .GT. 

GO TO 17. 

18 AMAX=D ( 1 » J ) 

NODE =J 
17 CONTINUE 

FORM THE FIRST SUBTOUR 

S< l,l)=NODE 
S!1,N0DE)=1 
LINKS! 1 )=2 

ESTABLISH ENDPOINTS 

53 CALL ENDPTS 

IF USER DESIRES TO PRINT SUCCESSOR FUNCTION AT EACH 
ITERATION SET IPRINT=3 

IFCIPRI NT .NE. 3) GO TO 54 
DO 109 1=1 » NTOURS 

WRITE(6,103) ( S( I, J) ,J=1,NC) 

108 FORMAT (2513) 

109 CONTINUE 

CALCULATE TOUR LENGTHS AND TOTAL DISTANCE 

54 CALL TORDIS 

ESTABLISH NONODE AND INNODE VECTORS 
CALL TOVIST 

WHEN I END=0 ALL CITIES HAVE BEEN VISITED 
I F ( IEND .EQ. 0) GO TO 50 
SELECT CRITERIA TO DETERMINE NEXT CITY TO VISIT 



999 STOP 

I F( NTOURS .GT. 
I F( NTOURS .EQ. 
CALL ENDPTS 
CALL THREEX 
GO TO 53 



1 ) 
1 ) 



CALL 

CALL 



GRATOR 

SNGLTR 



100 WRITE (6,105) ISTOP 

105 FORMAT ( 10X,'CITY* » I 5, ' IS OT OF RANGE*) 
GO TO 999 

50 DO 110 I = 1 , NTOURS 

WRITE(6,111) ( S( I, J) ,J=1,NC) 

111 FORMAT (2513) 

110 CONTINUE 

DO 112 1=1, NTOJRS 

WRITE! 6, 113 ) I , LENGTH ( I ) 

113 FORMAT! 10X,' LENGTH OF SUBTOUR 13 , • 

112 CONTINUE 

WRITE (6, 115) TOTDIS 

115 FORMAT! 10X, 'TOTAL DI STANCE F20 . 4 ) 

END 



= • , F 15 .4 ) 
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SUBROUTINE IN I T 



C SUBROUTINE WILL READ PROBLEM SIZE, CITY LOCATIONS AND 
C LOITER DISTANCES, AIRCRAFT RANGE, INITIALIZE SUCCESSOR 
C AND PREDECESSOR FUNCTIONS AND COMPUTE DISTANCE MATRIX 

COMMON D(25,25),S(25,25), INNODE ( 25, 2 5) , NONODE ( 25 ) , 
1RH(25 ,25) » LH( 25 , 25 ,L I NKS( 25 ) , THATA { 25 , 2 5 ) ,NTOURS,KK, 

II END, NC , TOTDi S, L ENGTH (25) , RANGE , I STOP , PRED ( 25 ,2 5) 
DIMENS ION X(25) .Y(25) ,STAR( 25,25) ,NNSTR( 25,25), 

ILL STR (25,25) , S=CSTR(25 ,25) , SNNS TR ( 25 , 25 ) , SL LSTRC 25 , 25 ) 
l,CRIT(25),OLDRH( 25,25) , OLDLH( 25,25) ,Z( 25) 

INTEGER S,RH,OLDRH,OLDLH, EXNODE, PRED 
REAL LENGTH 



I PRI NT=0 
I STOP=0 
NTOURS= 1 

C READ NUMBER OF CITIES IN PROBLEM AND AIRCRAFT RANGE 

READ (5,11) NC, RANGE 

11 FORMAT ( I3,F20. A) 

HALFR=RANGE/2. 

C READ CITY LOCATIONS AND LOITER DISTANCES, BASE MUST BE 
C READ AS CITY 1 

DO 20 1=1, NC 

RE AD ( 5 , 21 ) X( I),Y(I),Z(I) 

21 FORM AT ( 3F1 0* A’ ) 

20 CONTINUE 

C COMPUTE DISTANCE MATRIX 

C THE FACTOR 1.15 CONVERTS NAUTICAL MILES TO STATUTE MILES 

DO 12 I =1 , NC 
DO 13 J = 1 , NC 

D(I,J)=(SQRT(( ( X < I ) — X ( J ) )**2+{ (Y(I)-Y(J) )**2) ) ) )* 
1.15+Zt J) 

C CHECK FOR CITY OJT OF RANGE 

IF(I.EQ.l) GO TO 1 A 
GO TO 13 

1A IF ( D ( I , J ) .GT. HALFR) ISTOP = J 

13 CONTINUE 

12 CONTINUE 

C INITIALIZE SUCCESSOR AND PREDECESSOR FUNCTION 

DO 15 1=1, NC 
DO 16 J=1 , NC 
S( I, J)=0 
PRED ( I , J) =0 
16 CONTINUE 
15 CONTINUE 

C IF USER DESIRES TO PRINT DISTANCE MATRIX SET IPRINT=1 

IFUPRINT .NE. 1) GO TO 30 
DO 150 1=1, NC 

WR I T E ( 6 4 106 ) ( D ( I , J) , J = 1 ,NC ) 

106 FORMAT ( 25F A. 0) 

150 CONTINUE 
30 RETURN 
END 
C 

SUBROUTINE ENDPTS 
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on 



C SUBROUTINE ESTABLISHES, FOR EACH SU6T0JR, THE 

C ENDPOITS OF EACH LINK 



COMMON D< 25,25) ,S(25 ,25) . INNODE ( 25,2 5) , NONODE (25) , 
1RH(25 ,25) ,LH( 25 , 2 5 , L I NKS ( 25 ) , THATA ( 25 , 2 5 ) ,NTOURS,KK, 
1IEND,NC,T0TDIS,LENGTH(25) ,R ANGE, I STOP, PR ED ( 25,25) 
DIMENSION X ( 2 5) , Y ( 25 ) , STAR ( 25,25) ,NNSTR ( 25 , 25 ) , 
1LLSTR(25,25) , SECSTR( 25,25 ) , SNNSTR ( 25 , 25 ) , SL LSTR (25,25) 
I , CRI T( 25) , OLD RH( 25,25 ) , OLDLH( 25 , 25 ) , Z ( 2 5 I 
INTEGER S, RH, OLDRH, OLDLH, EXNODE, PRED 
REAL LENGTH 



I PRINT=0 

DO 22 J=I , NTOURS 
LH( J , 1 ) =1 
RH(J,1)=S( J,l) 

Li MI T=LI NKS ( J ) 

DO 23 K=2, LIMIT 
IM=K-1 
IL=RH( J , I M) 

LH ( J , K ) =RH( J » IM) 

IR=LH( J,K) 

RH ( J ,K) =S( J , IR) 

23 CONT IN J E 
22 CONTINUE 

C IF USER DESIRES TO PRINT THE ENDPOINTS OF EACH LINK AT 
C EVERY ITERATION SET IPRINT=4- 

I F( IPRI NT .NE. '+) GO TO 25 

DO IIO 1=1 , NTDJRS 
LIM IT=LINKS( I ) 

WRITE(6,102) (LHC I,J),RH(I,J) ,J = 1, LIMIT) 

102 F0RMAT(5X, 'ENDPOINT* ,2110) 

110 CONT IN JE 
25 RETURN 
END 

SUBROUTINE TORDIS 

SUBROUTINE CALCULATES LENGTH OF EACH SJBTOUR AND TOTAL 
DISTANCE AS SJM OF SUBTOJR LENGTHS 



COMMON D(25,25) , S ( 25 ,25 ), INNODE (25,25) , NONODE (25) , 
1RH(25 ,2 5) ,LH< 25 , 2 5 ,L I NKS ( 25 ) , THATA (25 ,25 ) , NTOURS, 
1IEND,NC ,TOTDIS, LENGTH (25) , RANGE, I STOP , 3 R ED ( 25 , 25 ) 
DIMENSION X(25) ,Y (25) , STAR (25,25) , NNSTR ( 25 , 25 ) , 

1LLSTR( 25,25) , S=CSTR( 25 ,25) , SNNSTR (25 ,25 ) , SL LSTR ( 2 5 , 25 ) 
1,CRIT (25) , OLDRH( 25 , 25 ) , OLDL H ( 25 , 2 5 ) , Z ( 2 5 ) 

INTEGER S,RH, OLDRH, OLDLH, EXNODE, PRED 
REAL LENGTH 



I PR I NT= 0 
TOTDIS=0.0 
DO 55 1=1, NTOURS 
NS=S( I, 1) 

NP=1 

LENGTH! I ) =0 . 0 

20 LENGTH! I )=D(NP,NS)+LENGTH( I ) 

NP = N S 
NS=S ( I , NP ) 

I F ( N P .GT. 1) GO TO 20 
TOTDIS=TOTDIS+LENGTH(I ) 

55 CONTINUE 

C IF USER DESIRES TO PRINT LENGTH OF EACH SUBTOUR AND 
C TOTAL DISTANCE AT EVERY ITERATION SET IPRINT=5 
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ODD 



IF ( IPRI NT .NE. 5) GO TO 25 
DO 105 1=1, NTOURS 

WRITE(6,102) I , LENGTH! I ) 

102 F0RMAT(13X, ' LENGTH OF SU BT3UR ' , 13 , ' =',F15.4) 

105 CONTINUE 

WRI TE ( 6 » 101 ) TOTDIS 

101 FORMAT! 10X, • TOTAL DISTANCE' , F20 .4 ) 

25 RETURN 
END 

SUBROUTINE TOVIST. 

SUBROUTINE ESTABLISHES A VECTOR (NONODE) OF THE CITIES 
THAT HAVE NOT BEEN VISITED AND FOR EACH SUBTOUR A VECTOR 
(INNODE) OF CITIES IN THAT SJBTOUR 



COMMON 0(25,2 5) , $ (25 ,25) , INNODE ( 25 , 25 ), NONODE (25) , 

1RH ( 25 , 25) ,LH( 25 , 2 5 , L I NKS( 25 ) , THATA ( 25 ,2 5 ) , NTOURS, 

II END, NC, TOTDIS, LENGTH (25) , RANGE, I STOP , PR ED ( 25 , 25 ) 
DIMENSION X ( 25) ,Y (25) ,STAR( 25,25 ) ,NNSTR( 25,25 ) , 

1LLSTR( 25,25) ,SECSTR( 25,25) , SNNS TR ( 25 , 25 ) , SL LSTR ( 25 , 25 ) 
1,CRIT(25) ,OLORH( 25,25 ) ,OLDLH( 25,25 ),Z( 25) 

INTEGER S, RH, OLDRH,OLDLH, EXNOD E,P RED 
REAL LENGTH 



KK=0 
I END=0 

DO 24 J=2 » NC 
DO 25 1=1, NTOURS 

IF ( S ( I , J ) .NE. 0) GO TO 24 

I F ( I .LT. NTOURS) GO TO 25 

IEND=1 

KK=KK+1 

NONODE ( KK) =J 

25 CONTINUE 
24 CONTINUE 

DO 26 1=1, NTOURS 
MM = 1 

DO 27 J=2 , NC 

IF ( S ( I , J ) . EQ . 0) GO TO 27 
I.NNO DE ( I , MM) =J 
MM=MM+1 
27 CONTINUE 

26 CONTINUE 
RETURN 
END 

SUBROUTINE SNGLTR 

C SUBROUTINE PLACES NEXT CITY IN SUBTOUR 



COMMON D( 25,2 5) , S ( 25 , 25 ) , I NNODE ( 25 ,25 ) , NONODE (25) , 
1RH(25 , 25) ,LH( 25,25,LINKS( 25) , THATA ( 25 ,2 5 ) , NTOURS, 

1 1 END, NC, TOTDIS, LENGTH (25) , RANGE, I STOP , P RED ( 25 ,25 ) 
DIMENS ION X( 25 ), Y( 25) , STAR ( 25,25) ,NNSTR ( 25,25), 

1L LSTR (25,25) , SECSTR ( 25 , 25 ) , SNNSTR ( 25 , 25 ) , SL LSTR ( 25,25) 
1 , CR I T ( 2 5 ) ,OLDRH( 25,25) , OLDLH( 25 , 25 ) , Z ( 25 ) 

INTEGER S»RH,OLDRH, OLDLH, EXNODE, PRED 
REAL LENGTH 



CALL TNCRMT 
AMAX=99999. 

AMI N=— 1.0 
I P ASS=0 

IF(LINKS( 1) . 3T . 4) GO TO 34 
LIMIT=LINKS(NTOJRS) 

28 DO 29 J=1 , KK 
DO 30 1=1, LIMIT 



ho 



IF (TH£TA( I » J ) .LT. AMAX) GO TO 31 
GO TO 30 

31 NN=MONODE(J) 

LL = I 

AMAX=THETA ( I t J) 

30 CONTINUE 

ST AR (1 » J ) = AMAX 
NNSTR ( 1,J) >NN 
LLSTR(1, J)=LL 
AMAX=99999. 

29 CONTINUE 

DO 32 K=1,KK 

IF(STAR(lt <) .GT. AMIN) GO TO 33 
GO TO 32 

33 IN=NNSTR( 1,0 
INL= LLSTR (1 , K) 

AM IN=STAR ( 1 » K ) 

32 CONTINUE 

IF C IPASS .EQ. 1) GO TO 35 
CHECK=L ENGTH( I ) +AMIN 
I F (CHEC K .GT. RANGE) CALL TOURNU 
I F ( CHEC K .GT. RANGE) GO TO 53 
I L=LH ( NTOURS » INL ) 

I R=RH ( NT OURS » INL) 

S ( NTOURS » I L ) =1 N 
S (NTOURS » I N ) = I R 
LINKStl ) = LINKS( 1)+1 
GO TO 53 



34 I PASS= l 
GO TO 28 

35 LIMIT=LINKS(NT0URS ) 

AMAX=99999. 

DO 36 I =1 » KK 
DO 37 J»l, LIMIT 

I F ( THETAt I , J) .GT.STARC 1, I ) .AND. THETA ( I , J ) . LT.AMAX ) 
1 GO TO 38 
GO TO 37 

38 AMAX = THETA( I , J ) 

LN K= J 

NOD=NONODE ( KK ) 

37 CONTINUE 

SECS TR ( 1 » I ) =AMAX 
SNNSTRt 1, I )=NOD 
SLLSTR ( 1 » I ) =LNK 
AMAX=99999. 

36 CONTINUE 

AMIN =- 1 . 

DO 39 K=1 . KK 

CR I T ( K ) =SEC S TR ( l f K ) -STAR ( 1 , K ) 

IF(CRITIK) .GT. AMIN) GO TO 40 
GO TO 39 

40 AMIN=CR IT ( K ) 

TEST =STAR( 1 » K ) 

NEWLNK=LLSTR(1,K) 

NEWNOD=NNSTR( 1»K) 

39 CONTINUE 

CHECK=L ENGTH ( 1 ) +TEST 
I F (CHECK .GT. RANGE) CALL TOURNU 
I F ( CHECK .GT. RANGE) GO TO 53 
I L = L H( NTOURS » NEWLNK) 

IR=RH( NTOURS tNEWLNK) 

S(NT0URS»IL)=NEWN0D 
S(NTOURS»NEWNOD)=IR 
LINKS (1 )=LINKS( l)+l 
53 RETURN 
END 

SUBROUTINE GRATOR 

C SUBROUTINE SELECTS NEXT CITY OR CITIES TO BE VISITED 
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C AND PLACES THE 'I IN THE APPROPRIATE SUBTOJR 



COMMON D( 25,2 5) , S(25 ,25) , I NNODE ( 25 , 25 ) , NONODE (25) , 
1RH(25,25) ,LH( 25 , 25 , L INKS( 25) , THATA ( 25 , 2 5 ) ,N TOURS, 

1 1 END »NC,TOTDIS,LENGTH(25) , RANGE, I STOP , PRED ( 25 , 25 ) 
DIMENSION X( 25) , Y( 25 ) , STAR ( 25,25) ,NNSTR ( 25 , 25 ) , 
LLLSTR(25,25) , SECSTR ( 25 , 25 ) , SNNSTR ( 25 , 25 ) , SL LSTR( 25, 25 ) 
1 , CR I T( 25 ) , OLDRH( 25,25) , OLDLH ( 25 , 25 ) , U 25 ) 

INTEGER S,RH,OLDRH, OLDLH, EXNODE, PRED 
REAL LENGTH 



I SAVTR=NT OURS 
DO 71 N =1 , ISAVTR 
NTOURS=N 
LIMIT=LINKS( N) 

CALL TNCRMT 
AMAX =99999 . 

DO 72 J=1,KK 
DO 73 1=1, LIMIT 

I F ( THETA( I , J ) .LT. AMAX) GO TO 74 
GO TO 73 

74 NN=NONODE(J) 

LL = I 

AMAX =THET A ( I , J ) 

73 CONTINUE 

STAR (N, J) =A MA X 
NNSTR ( N, J ) =NN 
LLSTR(N, J) =LL 
AMAX=99999 . 

72 CONTINUE 
71 CONTINUE 

NTOURS = ISAVTR 
I FLAG=0 

DO 75 I =1 , NTO JRS 
AMAX=99999 . 

DO 76 J= 1 , KK 

IF(STAR(I,J) .LT. AMAX) GO TO 77 
GO TO 76 

77 IN=NNSTR( I , J ) 

INL = LLSTR( I, J ) 

AMAX = STAR ( I, J) 

ICOMPR=J 

76 CONTINUE 

DO 78 K=1 , NTOURS 

IFU .EQ. I) GO TO 76 

I F ( AMAX .LT. STAR ( K , ICOMPR ) ) GO TO 78 
CHEC K=STAR ( < , I COMPR ) +LENGTH( K ) 

IF (C HECK .LT. RANGE) GO TO 75 

78 CONTINUE 

CHECK=AMAX+LENGTH( I ) 

IF (C HECK .ST. RANGE) 30 TO 79 
I FLA G=1 
I L=L H( I , I NL ) 

IR=RH( I , INL) 

S( I , IL)=IN 

S( I, IN ) =1 R 

LINKSt I )=LINKS( I ) +i 

GO TO 75 

79 STAR(I, ICOMPR) =99999 . 

75 CONTINUE 

I F ( IF LAG .EQ. 0) CALL TOJRNU 

CALL ENDPTS 

CALL TORDIS 

CALL TOVIST 

CALL CRSEXC 

RETURN 

END 

SUBROUTINE TNORMT 
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SUBROUTINE COMMUTES THE INCREMENTAL DISTANCE INCREASE 
OF PLACING AN* GIVEN NODE INTO ANY GIVEN LINK OF A 
SUBTOUR 



COMMON 0(25,25) , S ( 25 , 25 ) , I NNODE < 25 , 25 ) , NONODE (25) , 
1RH(25,25) »LH( 25, 25,LiNKS(25) , THATA(25 ,25 ) ,N TOURS, 

1 1 END ,NC,T0TDIS,LENGTH(25) , R ANGE , I STOP , P RED ( 25 , 25 ) 
DIMENSION X ( 2 5) , Y( 25 ) , STAR (25, 25) ,NNSTR ( 25 , 25 ) , 
ILLSTR(25 ,25) , SE C STR ( 25 , 25 ) , SNNSTR ( 25 , 25 ) , SL LSTR ( 2 5 , 25 ) 
1,CRIT(25) ,0LDRH(25,25) , OLDLH( 25 , 25 ) , Z ( 2 5 ) 

INTEGER S,RH, OLORH, OLDLH, EXNODE, PRED 
REAL LENGTH 



I PRINT = 0 

LIMIT=LINKS(NT3'UR$) 

DO 25 U = 1 * KK 
DO 27 1=1, LIMIT 
NL=LH( NTOURS , I ) 

NR=RH( NTOURS , I ) 

I N=NONODE ( J ) 

THETA( I ,U)=D(NL,IN)+D(IN,NR)-D(NL,NR) 

27 CONTINUE 
26 CONTINUE 

IF USER DESIRES TO PRINT COMPUTED VALUES OF THETA ( I , J ) 
SET I PR I N T=6 

I F ( IPRI NT .NE. 6) GO TO 50 
DO 112 1=1, LIMIT 

WRITE( 6, 113) ( THETA( I, J) ,J=1,KK) 

113 FORMAT (9F10. 2 ) 

112 CONTINUE 
50 RETURN 
END 

SUBROUTINE THREEX 



SUBROUTINE CHECKS CITIES IN A GIVEN SU3T0UR TO 
DETERMINE IF THE LENGTH OF THAT SUBTOUR MAY BE 
DECREASED BY RELOCATING THE CITY IN THAT SUBTOUR 



COMMON D( 25,2 5) ,S(25 ,25) , I NNODE ( 25 , 2 5 ) , NONODE (25) , 
1RH(25 ,25) ,LH( 25 , 25 , L INKSt 25 ) , THATA (2 5 , 2 5 ) , NTOURS, 
1IEND,NC,T0TDIS,LENGTH(25) .RANGE, I STOP , PR ED ( 25 , 2 5 ) 
DIMENS ION X ( 2 5) , Y ( 25 ) , STAR ( 25,25) , NNSTR ( 25 ,25 ) , 
1LLSTR(25,25) , SEC STR ( 25 , 25 ) , SNNSTR ( 25 , 25 ) , SLLSTR ( 25, 25 ) 
1 ,CRIT(25) ,0LDRH(25 ,25) , OLDLH ( 25 , 25 ) , Z ( 2 5 ) 

INTEGER S,RH,OLDRH, OLDLH, EXNODE, PRED 
REAL LENGTH 



55 DO 60 M=1 , NTO JRS 

IF( L INKS! M) .LT. 4) GO TO 48 
LIMIT=LINKS(M)-1 
DO 41 1=1, LIMIT 
NL=L H( M , I ) 

EXNODE=RH ( M , I ) 

NR=S(M, EXNODE) 

TSTVAL=D( NL, EXNODE )+D( EXNODE , NR) -D ( ML , NR ) 
INL IM= I +1 
NEWLI M=LI NKS ( M) 

DO L2 K= 1 , NE WL I M 

IF(K.EQ. I. OR. K.EQ.INLIM) GO TO 12 
I L=LH( M, K ) 

IR=RH(M,K) 

COMV AL =D ( IL, EXNODE )+D( EXNODE, I R ) -D ( IL , IR ) 
IF ( COMVAL .LT. TSTVAL) GO TO 45 
12 CONTINUE 



/ 




41 CONTINUE 

48 I F ( M .EQ. NTOURS) GO TO 46 
60 CONTINUE 
GO TO 46 

C CREATE THE IMPROVED TOUR 

45 S (M ,1 L) =EXNODE 

S ( M» EXNODE ) = IR 
S (M,NL) =NR 
CALL ENOPTS 
CALL TORDIS 
GO TO 55 

46 RETURN 
ENO 

SUBROUTINE TOJRNU 

C SUBROUTINE CREATES A NEW SUBTOUR WHEN IT IS NOT 

C FEASIBLE TO COMPLETE A VISIT TO EACH CITY WITH THE 

C PRESENT NUMBER OF SUBTOURS 



COMMON 0(25,25), S( 25, 25), INNODE ( 25, 25) , NONODE (25 ) , 

1RH (25,25) ,LH(25, 2 5,LINKS(25) , THA TA ( 25 , 25 ) , NT OURS, 

1 1 ENO »NC,TOTDIS, LENGTH (25) , R ANGE , I STOP , P RED ( 25 ,25 ) 
DIMENSION X ( 2 5 ) , Y ( 25 ) , ST AR ( 25 , 25 ) , NNS TR ( 25,25), 
ILLSTR(25 ,25) , SECSTR( 25 ,25 ) , SNNSTR ( 25 , 25 ) , SLLSTR( 25 , 25 ) 
L»CRIT(25)»QLDRH(25,25) , OLDLH ( 25 , 25 ) , Z ( 25 ) 

INTEGER S ,RH,OLDRH, OLDLH, EXNODc , PRED 
REAL LENGTH 



AMIN=99999. 

LIMIT=LINKS(N TOURS) 

NEWTDR=NT0URS+1 
DO 62 J =1 , KK 

NOD= NONODE ( J ) 

IF(D(I,NOD) .LT. AMIN) GO TO 63 
GO TO 62 

63 AM IN=D ( 1 , NOD ) 

I N=NOD 
62 CONTINUE 

S (NEWTOR, 1 ) = I N 

S(NEWTOR,IN)=l 

LINKS! NEWTOR) =2 

NTOURS=NEWTOR 

RETURN 

END 

SUBROUTINE CRSEXC 

C SUBROUTINE CHECKS CITIES IN A GIVEN SUBTOUR TO 

C DETERMINE IF TOTAL DISTANCE TRAVELED MAY BE DECREASED 

C BY PLACING THEM IN A DIFFERENT SUBTOUR 



COMMON 0(25,25) ,S (25 ,25) , INNODE ( 25, 25 ), NONODE (25 ) , 
IRH(25,25) ,LH(25,25,LINKS(25) , THATA ( 25 ,25 ) ,NTOURS, 

1 1 END, NC, TOT D IS, LENGTH (25) , RANGE, I STOP , P RED ( 25 , 25) 
DIMENSION X ( 2 5 J , Y ( 25 ) ,STAR(25,25) . NNSTR ( 25 , 25 ) , 
ILLSTR(25,25) , SECSTRI 25 , 25 ) , SNNSTR ( 25 , 25 ) , SL LSTR ( 2 5 , 25 ) 
I ,CRIT( 25 ) ,0LDRH(25,25) , OLDLH (25 ,25 ) , Z( 25 ) 

INTEGER S,RH,OLDRH, OLDLH, EXNODE, PRED 
REAL LENGTH 

86 DO 80 1=1, NTOURS 
NOD=S ( I ,1) 

PRED( I , NOD) =1 
81 INT=S( I ,NOD) 

PRED(I , INT ) =NOD 
NOD= I NT 



hh 



IF ( I NT .NE. 1) GO TO 81 
80 CONTINUE 

DO 82 I=l,NTOJRS 
LIMI T=LINKS ( I ) 

DO 83 J=1 » NTO JRS 

IF C I .EQ. J) GO TO 83 
N JMB ER = LI N<S ( J )-l 
I F ( NUMBER .EQ. 1) GO TO 82 
DO 84 K=l, NUMBER 
DO 85 M=l, LIMIT 
NL=LH( I t M) 

NR=RH( I ,M) 

IN= I NNODE ( J t K ) 

IBFR=PR£D! J, IN) 

IAFT = S( J, I VI ) 

THETA! M,K) =}( NL, IN )+D( I N , NR ) -D (NL , NR ) 
BETA=D( IBFR, I N ) +D ( IN , I AFT ) -D ( I B FR, UFT) 
I F ( T HET A( M » < ) .GE. BETA) GO TO 85 
CHECK=THETA( M,K)+LENGTH( I ) 

IF ( C HECK . 3T . RANGE) GO TO 85 



S< J, IBFR-) = IAFT 
S( J, IN) =0 

LINKS! J)=LINKS(J)-1 

INSERT IN TOUR I 

S( I » NL) = 1 N 
S( I , IN ) =NR 
LINKS! I ) = L INKS ( I ) +1 
CALL ENDPTS 
CALL TOVIST 
CALL TORDIS 
GO TO 86 
85 CONTINUE 
84 CONTINUE 
83 CONTINUE 
82 CONTINUE 
RETURN 
END 
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